Gray soliton solution in the extended nonlinear Schrodinger equation 
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In the framework of extended nonlinear Schrodinger equation (ENSE) the classification of self- 
similar solutions by the relation between the amplitude and phase is performed. New solutions of 
ENSE - "gray soliton" and "gray soliton chain" are presented. The properties of these solutions and 
the possibility for using theirs in physical applications are discussed. 

PACS numbers: 42.65.-k, 75.30.D 
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The evolution of the waves in a weakly nonlinear higher 
order dispersion media can be described by means of the 
extended nonlinear Schrodinger equation (ENSE) 
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The first three terms in form "classic" nonlinear 
Schrodinger equation (NSE). The coefficients of equation 
are connected with properties of propagating waves. So, 
a ~ uj'lf. is the group velocity dispersion (w is frequency 
and k is wave number of carrier wave), 0:3 describes 
the third-order dispersion, N connects with nonlinear re- 
sponse of the medium. The coefficients ai, a2 describe 
the nonlinear dispersion characteristics of medium and 
their nature depends on the problem under consideration. 
For example, the appearance of such terms in nonlinear 
optics is caused on inhomogeneous raman scattering (see, 
for example Q), in the nonlinear ferromagnetodynamics 
the values ai = a2 show the dependence of nonlinear 
response of medium from wave number. 

The Eq. is not completely integrable in terms of 
inverse scattering problem, but there are some integrable 
cases. The simplest of them ai = a2 = az = 0, when 
equation reduce to NSE. The another one is ai = 
a2 = 6a3, in this case reduces to modified Korteweg 
- de Vries equation by means of Galileo transformation. 
For Q!2 = and aai — Sa^N (Hirota conditions), the 
ENSE reduces to the Hirota equation 1^ and for a — 1, 
N = 2, 2ai = 80:3, 4ai = 803 eq. ^ reduces to the 
Sasa - Satsuma equation. 

In general case (with arbitrary values of parameters) 
equation 0J is not completely integrable, but some solu- 
tions (and even soliton-like) are known. In our opinion, 
the most important of them are the follows: "light" and 
"dark" Potasek - Tabor solitonsja] , , so-called "embed- 
ded" and "radiating" solitons 0, cnoidal states Q and 
some specific solutions (including "soliton with tail" and 



"algebraic soliton" ) which were found in ^ (but these so- 
lutions realize only for specific values of parameters ai. 

In this paper we try to extend the class of exact so- 
lution of ENSE in non integrable case (with arbitrary 
values of coefficients) . We show here that a wide class of 
solutions can be classified in form of the relation between 
the nonlinear frequency shift and amplitude of solution. 

We suppose solution of (Q) in the form of a stationary 
wave with phase modulation 

fiv, t) = F{y - vt) exp{i(TO ~qt + a{y - vt))}. (2) 

In this case the equation (QJ reduces to the system of two 
ordinary differential equations (ODEs) for real functions 
F{x), Sl(a;) (we denote here x — y — vt, n{x) — a'^) 



(a - 3a3[il + p))—r- - 80:3 — — 
dx"^ dx dx 

+ [N - ai{p + n)]F^ + 



(3) 



d^n 



+ [-a{n+pY + asin +py-a3-^ + nv + q]F = 

d^F dF 
as-^ + [2a(n+p)-3as{n+p)'-v]—+ (4) 



(ai + 2a2)F^^ + (a - Sa^ip + fl))F^ 
dx dx 



0. 



Let us consider the solution of the system of equations 
©, 101) in the following form 



n = A + BF" 



(5) 



Note, that we do not impose restrictions on possible val- 
ues of n - it may be any real number, both positive and 
negative, integer and fractional. Note also, that case 
n = 1 was discussed in 9] and the case B = was consid- 
ered in 01 . By substituting (|SJ| into the system of ODEs 
(O, 10} and integrating the equation 0} with respect to 
X we come to the system 
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(a -3a3{A+p))-j-j + [N - ai (A + p)]F^ + [-a{A + pf + aaC^ + pf + vA^ q]F 



(6) 



H^F / 

in + 3)a3BF''^ - n(n + 2)a3SF"-i — 
ax \ ax 



[-2a(A+p) + 3a3{A + pf + v]BF^+'' - aiBF^+'' ■ 



\l+3n 



+ i-a + SaaiA + p))B' F{xy+''^ + a^B^Fix) 

d?F 1 
= as-ij^ + [2a{A + p)-iaz{A+pf ~v]F + -{ai+2a2)F 

1 + 2n 



n + 2 
n+1 



{a-3a3{A+p))BF 



l+n 



(7) 



Let's compare the Eq. ^ and Eq. {Tj). The first of them 
has two "bad" terms: 

-{n + 3)a3BF"F^^, (8) 

which is proportional to the product of n power of func- 
tion F by its second derivation, and 

-n(n + 2)a3i3F"-i(F^)2, (9) 



containing the product of rt — 1 power of F and square 
of its first derivation. We can easy get rid these terms in 
the following way. At first we multiply the equation (UJ 
by {n + 3)BF" and add the result to the ©■ Second, we 
multiply Eq. (O by F^, integrate the result, next find 
the relation for (F^)'^ from the obtained equation and 
substitute the result into ©. Thereafter we have the 
system of two ODEs of the same form: 



[a ~ 3a3{A + p))^ + [~a{A + pf + a^iA+pf ^vA + q]F + [N - ai{A+p)]F^ + 2n{n + 2)C2BF^-^ 
{n^ + 3n + 2)[2a{A + p) - 3a3{A + pf - v]BF^+'' + {3n + 5){a - 3a3{A + p))B^F^+^'' + 
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(2n" + 5n + 3)SCiF" + [(ai + 2a2)— + (ai + 2a2)— + 2a2]BF^+'' - 2 

6 3 

1 n + 2 

as-TY + + P)- ^MA+pf - v]F + -(ai + 202)^^ + — - 
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d^F 

' dx'^ 
l + n 

l + 2n 



3n2 + 8n + 4 
1 + 2n 



asB-'F 



3 T7'l+3n 



{a-3a3{A+p))BF 



l+n 



a^B^F^+^" +Cu 
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(10) 
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where Ci and C2 are constants of integration. The sys- 
tem of this ODEs has nontrivial solutions if the equations 
(|10|l . coincide, that is they have to contain the value 
F{x) with same indices of power and proportional co- 
efficients. In this context there are two cases available: 
(i) all exponents in pU|l . are different and (ii) some 
values coincide. 

(i). If all exponents have different values, the condi- 
tions of compatibility of the system 110() . have the 
follows form: 



: Ci = 

^"-1 : C2 = 



(12) 



F 

p3+n . ^ 

: aaS^pn^ + 8n + 4] = 



TL ^Jl 

— {ai + 2a2) + -;r("i 
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We don't consider here the case _B = - in this case 



we immediately come to the Duffing equation and brief 
investigation of available solutions in this situation can 
be found in 8] . It is clear from last two equations in 112|) 
that system with B ^ Q may have nontrivial solutions in 
the only one case: 



ai 



-2/3, 
17/5a2. 



(13) 
(14) 



So we have remarkable result: extended NSE has non- 
trivial solutions only for bounded set of numbers n. The 
equation for n = —2/3 has the follow form: 



= as 



d^F 
dx"^ 



9 



a2F^ + 9(a - 3a3p)fcoi^ 



(15) 



4(a - 3a3p)BF^/^ + 3a3B^F-^ 



/3 



3 



with parameters 

WasN ~ 17aa2 



U(F) 



kn 



A + p 



10q;2Q^3 



(16) 



10a2a3 

V = 2a{A+p)~3a3{A+p)'^ -9{a-3a3p)ko, 
q = a{A+pf -a^iA+pf - Av + 
+ ko{2a{A + p) -3a3{A + pf -v). 

At first sight the existence of physical system with singled 
out parameters (|14|1 is low - probability. 

(ii). Let's investigate now the case, when there are 
some coincident values among the indices of power of 
function F{x) in H1U|) . Hll|) . It appears when n takes on 
one of the values 



0, 1, 2, 3, 4, -1, 
-1/2, -1/3, 2/3. 



(17) 



For each value from H17|) we have its own system of equa- 
tions and all cases should be analyzed separately. Note, 
the case n = is identical to the case B = (we can re- 
name A + B Ain such case). The values n = —1,-1/2 
are specific because of zeros in denominator. For this val- 
ues we are to do full transformation chain Q - (|ll|l . As 
for other values from H17|) the final system of ODEs can 
be found after direct substitution the corresponding value 
of n into (UHl) and ifTTl) . 

After do this we come to the follow conclusion. Non- 
trivial solutions may be found only in three cases: 

1. B = 0. Duffing equation. The classification of pos- 
sible solutions was performed in 

2. n = 1. The brief classification of possible states in 
this situation was published in 9] . 

3. n = -2. This case will be discussed bellow. 

The substitution n = 2,3,4,-1,-3,-1/2,-1/3,2/3 re- 
quires the condition = 0. As result we have instead 
of third-order differential equation a second-order ODE. 
We are not interesting this situation in our paper. 
In the case n = — 2 we have the equation 



dx 



2 , a^^{2a{A^p)-?,a:i{A + pf -v)F + (18) 
+ ll3a-^{ai + 2a2)F^ - B^F-^ = 



whose parameters satisfy the relations 

Ci = 0, 

a - 3q3(A-|-p) _ N ~ai{A+p) 

— o ~ 



(19) 



{a2 - ai)B - a{A + p)'^ + asiA + p)^ + vA + q 
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FIG. 1: Phase portrait for the case (ai + 202)03 < 



The first integral of motion for this equation can be easy 
calculated 



-aj\2aiA+p) - 3a3{A + pf - v)F'^ - 
l/4,a^^{ai+2a2)F'^ - B^F-^ + 2E. (20) 



Let's consider now the case [ai + 2^2)0.3 < 0. In this 
situation potential (right - hand part of H2()(l l is negative 
for small values of F and increases infinitely for large 
F. So, for small enough F equation H18|l have no real 
solution, whereas for large F the solutions have infinite 
trajectories. For some values of E potential has three real 
roots and the finite periodical solutions are possible. This 
potential for different values of E and corresponded phase 
portrait is shown in fig. ^ The solutions, corresponded 
to the finite trajectories in phase portrait on fig. ^ can 
be defined by means of Jacobi elliptic delta function 



2a{A+p) - 3q!3(^ +p)2 - V 



2 

Fix) = Fo[l-^^dn{bx,kf} 



1/2 



(21) 
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FIG. 2: "Gray" soliton in ENSE model 



with parameters 



Fo 



E = 



2?7(l + /i2) ^1/2 



3 + /i^ 



-,1/2 

2^2 + 2fcV^ + 3 



(22) 



^2 (3 + ^2_^fc2^2)3 



where ^, 77 are coefficients standing before F'^ and F in 
(|18|l respectively. The separatrix solution (potential has 
the point of contact with F axis, see fig. ^ appears when 
fc = l(0<fc<l - parameter of elliptic Jacobi function) . 
Returning to the space variable y and to the function 
ip{y,t) and using the relations (|22|l . p9|) we find gray 
soliton solution of extended NSE 



F{y~vt) = 
cr{y - vt) = 



F{y - vt) cxp{i{py 
,2 



^^0 1 



-sech" 



f2 



1 + 1.1 

arctan[/i tanh(6(2; — vt))] 
6a3 1 



-qt + a{y -vt))}, 
1 1/2 

(Hy-vt))} , 



^%2^ 



q;i + 2a2 A^2 
a(ai + 2a2) — SNa^ b 
60302 A*' 



(23) 



= ap2 _ Q-gp^ - F^{N - Qip), 

= 2q!P — 3ck3p^ — 2036^ — 2{a — 3a3p)b/iJ. 

— Qa^iP' / . 



The gray soliton solution of "classic" nonlinear 
Schrodinger equation was known long ago (see for in- 
stance [13) but such solution in the ENSE was not 
known. The evolution of this solution is shown in fig. 
121 For numerical simulation we use Fourier transform 
in y-space and fourth-order Runge-Kutta method in t- 
space. 

In spite of Potasek - Tabor solitons |^ this solution 
has two free parameters (for example thickness and 
deep of modulation /{\ + y?), < n < 00). Note 
also, that amplitude and phase in this solution are con- 
nected by means the relation So we can control the 
deep of modulation by means of nonlinear phase shift, 
and otherwise we can generate different phase shift using 
the corresponded deep of modulation. These properties 
should be fruitful for different physical application. 

As for cnoidal solutions 1)21(1 at A; 1 they represent 
well separated holes on against a background of the car- 
rier wave and can be understand as "ID gray soliton lat- 
tice". Here we discuss the case {ai + 2a2)ct3 < 0. In the 
opposite case equation H18|) also has finite periodical so- 
lutions but this situation should be discussed separately. 
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